Effect of the surface polarization in polar perovskites studied from first principles 
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The (001) surfaces of polar perovskites BaTiOa and PbTi03 have been studied from first principles 
at T=0 K. For both cases of polarization, the most stable TiO-terminated interfaces show intrinsic 
ferroelectricity. In the topmost layer, where the O atoms are >0.1 A above Ti, this leads to metallic 
instead of the insulating behavior of the electronic states that may have important implications for 
multiferroic tunneling junctions. 
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Epitaxial growth technique opens the way to combine 
ferroelectrics (FE) and ferromagnets (FM) into multi- 
component multiferroics, which are inaccessible by tra- 
ditional synthesis* 1 -! 2 - Currently, many laboratories across 
the world work on a prototypic device, consisting of a 
few-nm thick ferroelectric sandwiched between a ferro- 
magnet and another metallic contact. External electric 
field applied to the FE phase causes a switching of its 
polarization P that, in turn, through the magnetoelec- 
tric coupling may change the magnetic order in the fer- 
romagnetic phase*? Altering P and magnetization inde- 
pendently for encoding information in multiferroics, the 
smallest quaternary logical memory might be obtained. 
Fundamentally new multiferroics require a better under- 
standing of the FE order parameters, especially, in the 
case of thin films here the symmetry is reduced. 

Since direct measurements of atomic displacements, 
occurring in FE near the interface, is extremely chal- 
lenging, their structures can be understood and nu- 
merically characterized from first principles. Recently, 
much work has been conducted to study bare FE sur- 
faces using the ab initio density functional theory (DFT) 
calculations4i5iSiLSi2ii2iii i i£ It has been found that the 
critical thickness down to 3 unit cells (1.2 nanometers) is 
enough to enable the existence of ferroelectricity at room 
temperaturei 13 i 14 ' 15 However, there are only some con- 
vincing evidences in literature i.e. the work by Cohen 5 
that the direction of P may affect the surface relaxation. 
The functionality of multiferroics assumes that P must 
be reversible and parallel to the surface normal. Hence, 
it is worthwhile to carry out ab initio calculations which 
model the reaction of the (001) surface of polar ferroelec- 
tric surfaces upon the change of its reversible P. 

In this report, we study the (001) surface of ABO3 
perovskites (A = Sr, Ba, Pb and B = Ti), which repre- 
sent a wide class of ferroelectrics ranging from paraelec- 
tric SrTi0 3 (STO) to highly polar PbTi0 3 (PTO), while 
BaTiOa (BTO), with its moderate spontaneous polar- 
ization P s , is an example of a typical FE. The study is 
based on extensive calculations, using the Vienna Ab- 
initio Simulation Package 16 (VASP), in which the effects 
of relaxation of atomic positions are included. Nowadays, 
many FE properties can be successfully calculated from 



first principles.— Table Q] collects the experimental data 
for the lattice parameters and atomic positions, obtained 
for the room-temperature tetragonal phase of PTO and 
BTO (with space group symmetry PAmm) and for cu- 
bic STO (Pm — 3m), in comparison with our calcula- 
tion. Overall, there is a good agreement between the 
measured and theoretical structure parameters for the 
three systems. For polar BTO and PTO, their values of 
P, calculated by the Berry phase approac h 22 ' 23 , are in 
reasonably good agreement with experiment. The minor 
differences, seen in Table U between our and some other 
recent DFT results* 6 - ' 12 ' 24 can be attributed to the choice 
of pseudopotentials and/or to the approximation used 
for the exchange and correlation potential. We used the 
local density approximation (LDA) while the electron- 
ion interactions were described by the PAW pseudopo- 
tentials. After relaxation the calculated forces are al- 
ways less than 0.5-10 _2 eV/A. The electron pseudo-wave- 
functions were represented using plane waves, with a cut- 
off energy of 650 eV. For the Brillouin-zone integration a 
dense Monkhorst-Pack^ 5 - mesh was used. 

We calculated 5-unit-cell (~2-nm) thick ABO3 film. 
The atoms of the two upper or, alternatively, two lower 



TABLE I: The lattice parameters a, c/a, fractional atomic 
coordinates, z/a, and polarization, calculated for the room 
temperature PTO, BTO and STO phases, are shown in com- 
parison with the corresponding experimental data. 
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unit cells were allowed to relax while all other atoms in 
the supercell were fixed at their bulk-like and previously 
optimized positions, which are shown in Table HI A vac- 
uum spacer of 2 nm was used to separate the copies of 
the periodic structures in the direction perpendicular to 
the surface. The ABO3 perovskite structures possess 
a strong anisotropy resulting in the AO (A = Pb, Ba) 
and TiC>2 layers alternating in the [001] direction. The 
(001) surface of a ABO3 can be terminated by a AO or 
a Ti02 layer. Recently, Eglitis and Vanderbili— have re- 
ported for the cubic structure of BTO and PTO that the 
Ti02-terminated surface is more stable. Using the same 
approach^ we calculated the surface energy for the both 
terminations of ABO3. The results are shown in Table HT1 
For each perovskite, its Ti02 terminated (001) surface is 
energetically favourable, indeed. Thus, we consider the 
(001) surface of ABO3 to be Ti02 terminated in the fol- 
lowing. 

In polar PTO and BTO, the displacements of the Ti 
and O atoms occur along the z-axis so that P is consid- 
ered to be directed along [001] direction as well. First, 
we must formally set up the direction of electric dipole 
in the unrelaxed supercell assuming, for instance, that O 
is always above the corresponding cations in each layer 
along [001], as given in Table fl] Then we can model the 
two distinct situation alternatively placing the bulk-like 
3-unit-cell thick substrate and relaxed layers against each 
other along the z axis. In the first case, which we denote 
as Pi, the direction of P is antiparallel to the surface 
normal. The second model labelled by Pf corresponds to 
the case where all cations are above O before relaxation 
and, therefore, P is parallel to [001]. In the tetragonal 
FE structure, both configurations may coexist in the ran- 
dom state as Pf and Pf domains separated by a domain 
wall of <2 nm. To quantify the process of relaxation we 
use the cation-anion displacements 5 = zo — z cat i on cal- 
culated for each AO and TiO layer near the interface. 
In the bulk- like substrate of polar PTO and BTO, the 
model Pf means that 8 < and vice versa the case of Pi 
models the situation where 6 > 0. 

Fig. rjjshows several top monolayers (ML) of PTO(001) 
after relaxation. The case Pf (Pf) is shown on the left 
(right) side of Fig[TJ The arrows indicate the direction 
of the dipoles in each ML, while the numbers at the ar- 
rows give the intralayer displacements 5 in A , calculated 
between the O and metal atoms along [001]. In the state 
Pi, all dipoles possess the same orientation that means 
that O is always above the cation within each layer. In 
bulk PTO, the intralayer displacement, 5, is 0.333A for 
the TiO layer and 0.476A for the PbO layer. For the 
three top ML near the interface, their 6 are reduced by 
30-45% with respect to the corresponding bulk values. In 
the topmost TiO layer, the reduction of S is ^30%. For 
the second (PbO) and third (TiO) layers from the inter- 
face, we find that their 5 is reduced by 45% and 40%, 
respectively. In the case of Pf , shown in the left panel of 
Fig. [U the result of relaxation is rather different. In the 
third ML, the separation between Ti and O along [001] 
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FIG. 1: The relaxed top layers of the PTO (001) surface, with 
polarization Pf (left) and Pf (right). The arrows indicate the 
direction of the dipole within each layer while the numbers 
show the displacements 5 in A between O and Ti/Pb along 
[001]. For the case of Pf the sign of 5 in the topmost TiO 
layer is reversed. 



is 0.153 A , which is reduced by 54% against the corre- 
sponding bulk value. For the second ML, we obtain the 
reduction of 33%. However, the most significant changes 
occur in the topmost TiO layer, whose S is largely re- 
duced by 68% whereas the dipole is reversed compared 
to all others. Thus, using the Pf model and placing all 
O below the cations, we obtained in the topmost layer 
the relaxed configuration where O is above Ti. This is 
similar to the case of Pf . 

To investigate the effect of the surface rumpling in 
perovskites we repeat the calculations for STO(001) and 
BTO(001). In Table HT1 the corresponding results of our 
zero temperature calculation are listed. For paraelectric 
STO, we obtain that its TiO-terminated (001) surface 
after relaxation becomes marginally polar, with a posi- 
tive rumpling normal to the surface in the top three ML 
where the O atoms are above the cations by <0.12 A . 
This is in good agreement with the most recent exper- 
imental studies^ 6 The positive rumpling predicted for 
bare surfaces of perovskites leads to relatively low cat- 
alytic activity. With increasing temperature, the rum- 
pling is distorted and it may stimulate further potential 
catalysis. For the TiO-terminated BTO surface, we have 
found the details of relaxation similar to those of PTO. 
In fact, all our results are in good agreement with those 
reported by Eglitis and VanderbiltJ^ In case of Pf, the 
topmost BTO rumpling of ~0.1 A being larger than the 
corresponding bulk value, is similar to that of highly po- 
lar PTO. The sign of S in the topmost ML of BTO is 
reversed with respect to all others calculated for the lay- 
ers situated far down from the interface. In the third 
ML S is 0. The Pf model yields for BTO(001) the re- 
versal dipole in the second ML, with marginal 5. Thus, 
we find for the three systems and different arrangements 
of P that the TiO-terminated (001) surfaces prefer the 
configuration where O is above Ti. In the cubic ABO3 
perovskite structure, each Ti 4+ ion sits in the regular six- 
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TABLE II: The surface energy E surf (in eV) for the Ti0 2 
and AO terminated (001) surfaces and the cation-anion dis- 
placements 5 (in A) calculated for the top five layers of the 
perovskite (001) surfaces. 
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FIG. 2: Ti-0 bond lengths in the fivefold coordinated poly- 
hedra of the topmost layer of PTO(001), BTO(001) and 
STO (001) compared to those obtained in the corresponding 
bulk structures, with the sixfold coordinated environments. 



fold coordinated site with all of the Ti-0 bonds of equal 
length, as shown in Fig. [2] for bulk STO. 

In the tetragonal perovskite structure, such as t-PTO, 
the relaxed cluster of O atoms about the sixfold coordi- 
nated Ti forms a distorted octahedron, where one of the 
two bond lengths along [001] is rather short while an- 
other Ti-0 bond in the vertical direction is significantly 
longer than the four other bonds stretching in the equato- 
rial plane. If we exclude the longest Ti-0 bond from the 
consideration using the electrostatic arguments then the 
environment for each Ti becomes the fivefold coordinated 
polyhedra, which is similar to that of Ti at the interface. 
The left three panels of Fig. [2] compare the Ti-0 bond 
lengths in bulk PTO, normalized to the value of the ideal 
octahedron in the cubic structure, to those in the polyhe- 
dron around fivefold coordinated Ti in the topmost layer. 
Using the Pj model for the PTO (001) surface, we ob- 
tain the Ti-0 bonds whose lengths are similar to those of 
i-PTO and, hence, nothing dramatic happens in the en- 
vironment of the topmost Ti. In the case of Pj, the bond 
length distribution around the topmost Ti is restricted so 
that the equatorial and vertical bond lengths tend to be 



equal to each other. Moreover, the closest O atom to the 
surface Ti is attached along [001] from the opposite side 
compared to all Ti placed below the surface in the regular 
crystal structure within Pj. It is clear that the O-Ti-0 
bond angles for the equatorial Ti-0 bonds of the top- 
most ML must be dramatically changed to compensate 
the charge distribution around Ti. It appears that these 
bond angles become >90°, as shown in Fig. [TJ There- 
fore, whatever the state of P is modeled in i-PTO, the O 
atoms must relax above Ti on the Ti-0 terminated (001) 
surface. Regarding BTO, the same conclusions may be 
drawn. 

Fig. [3] shows the view of the charge density along [010] 
calculated for the top six ML of PTO and projected on 
the x-z plane of the supercell. The isocharge lines plotted 
in Fig. [3]for both cases of P illustrate the charge transfer- 
across the cell while the arrows indicate the dipole direc- 
tions within each ML. In the case of Pj, which is shown 
in the right panel of Fig. [3j there are three bridges seen 
between Ti and nearest O. The shortest bond with O, 
which is always below Ti along [001], has the large popu- 
lation value. In the Pj state, the charge transfer picture 
is similar to that of Pj for the topmost Ti only. Far be- 
low the interface (starting from the 5th ML) all Ti have 
the most populated bond with O which is above Ti. In 
the third ML, however, the Ti ion is strongly bonded to 
equatorial oxygens showing some sort of blockade for the 
charge transfer along the [001] direction. This may re- 
veal the key electronic states factors behind the surface 
relaxation of polar FE. 

Recently, Urakami et al.— have observed the surface 
conductance of BaTiOa single crystals in ultra high vac- 
uum below Tc- It has been shown that the in-plane con- 
ductance is the result of an intrinsic surface electron/hole 
layer that is, due to the surface polarity and not due to O 
vacancies or some other defects. The I — V characteris- 
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FIG. 3: The charge density, calculated for the top six layers 
of PbTiO3(001) and projected on the x-z plane. The two op- 
posite polarization states are shown: Pj (left) and Pj (right) 
where the dipole directions within each layer is labelled by 
the arrows. In the case of Pj , the isocharge lines of the third 
TiO layer from the interface show a blockade of the charge 
transfer along [001]. 
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FIG. 4: The Ti- and Oi-resolved DOS of t-BTO are plotted 
in the two upper panels where the top of the valence band is 
taken as zero energy. In the lower panels, the corresponding 
local DOS of the TiO-terminated (001) surface are plotted for 
the Pf and Pf states. 



approximation of DFT. Comparing the Ti and O DOS 
of t-BTO and the topmost ML of BTO(OOl), we see a 
spectacular change of the electronic states occurring due 
to the surface relaxation and variation of P. The major 
DOS features can be summarized as follows. In the case 
of Pf , a few O states appear in the band gap while the Ti 
DOS is not affected. For the Pf poled state, the Ti lower 
conduction band, being shifted downwards in energy by 
~2 eV, significantly contributes to the DOS in the band 
gap region. This causes metallic behavior of the topmost 
ML in the case of Pf , yielding rather large in-plane con- 
ductance. In the case of Pf the Ti states have a gap at Ej 
which is related to a tiny in-plane conductance. Depend- 
ing on the polarization direction the topmost ML under- 
goes a transition from metallic to oxide behavior shows 
metallic or oxide behavior. As a consequence the in-plane 
conductance changes drastically which is a reasonable ex- 
planation of the experimental results by Urakami et al.— . 



tics shows a pronounced difference of conduction between 
the poled states in BTO. We can explain this difference 
in a simple way using our ab initio results. To reveal 
the differences between Pf and Pf we plot in Fig. 2] the 
Mulliken site-projected density of states (DOS) of the 
BTO(OOl) surface for both cases of polarization. The 
Ti and O DOS for the topmost ML are shown in com- 
parison with the corresponding DOS of i-BTO. For bulk 
BTO, Fig. 0] shows a pronounced insulating band gap of 
2 eV. The value is typically underestimated by the LDA 



In summary, from the ab initio basis of our work we 
have shown that the intrinsic ferroelectricity in polar per- 
ovskites is suppressed by ~30% in the surface region. 
For both cases of polarization direction, the TiO termi- 
nated surface of BTO and PTO forms an electric dipolc 
where the O atoms being shifted >0.1 A above Ti. But 
nevertheless the electronic structure of the surface layer 
changes from metallic to oxide behavior under reversal 
of polarization which changes the surface conductance 
drastically^ This may have important implications in 
the design of multiferroic nano-devices. 
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